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ABSTRACT 

Context. The INTEGRAL mission has led to the discovery of a new type of supergiant X-ray binaries (SGXBs), whose physical properties 
^> . differ from those of previously known SGXBs. Those sources are in the course of being unveiled by means of multi-wavelength X-rays, 
f*") optical, near- and mid-infrared observations, and two classes are appearing. The first class consists of obscured persistent SGXBs and the 

second is populated by the so-called supergiant fast X-ray transients (SFXTs). 
f" — ■ ■ Aims. We report here mid-infrared (MIR) observations of the companion stars of twelve SGXBs from these two classes in order to assess the 
H contribution of the star and the material enshrouding the system to the total emission. 

, Methods. We used data from observations we carried out at ESO/VLT with VISIR, as well as archival and published data, to perform 
f""*) ' broad-band spectral energy distributions of the companion stars and fitted them with a combination of two black bodies representing the star 
OO an d a MIR excess due to the absorbing material enshrouding the star, if there was any. 

O ■ Results. We detect a MIR excess in the emission of IGR J16318-4848, IGR 116358-4726, and perhaps IGR J16195-4945. The other sources 
do not exhibit any MIR excess even when the intrinsic absorption is very high. Indeed, the stellar winds of supergiant stars are not suitable 
. _ for dust production, and we show that this behaviour is not changed by the presence of the compact object. Concerning IGR 116318-4848 
■ and probably IGR 116358-4726, the MIR excess can be explained by their sgB[e] nature and the presence of an equatorial disk around the 
?H ' supergiant companion in which dust can be produced. Moreover, our results suggest that some of the supergiant stars in those systems could 
5^ i exhibit an absorption excess compared to isolated supergiant stars, this excess being possibly partly due to the photoionisation of their stellar 
wind in the vicinity of their atmosphere, (abridged) 



Key words. Infrared: stars - X-rays: binaries, individuals: IGR 116195-4945, IGR 116207-5129, IGR J16318-4848, IGR J16320-4751, 
IGR 116358-4726, IGR 116418-4532, IGR 116465-4507, IGR J16479-4514, IGR J17252-3616, IGR 117391-3021, IGR J17544-2619, 
IGR J19140+0951 - binaries: general - supergiants - Stars: fundamental parameters 



1. Introduction 

High-mass X-ray binaries (HMXBs) are X-ray sources for 
which high-energy emission stems from accretion onto a com- 
pact object (black hole or neutron star) of material coming from 
a massive companion star. Until recently, the huge majority of 
known HMXBs were Be/X-ray binaries, i.e a neutron star ac- 
creting from a disc around a Be star. Most of these sources 
are transient, even if a few are persistent weak X-ray emitters 



Send offprint requests to: F. Rahoui 

* Based on observations carried out at the European Southern 
Observatory under programmes ID 075.D-0773 and 077.D-0721. 



(L x ~ 10 erg s ). The other known HMXBs were super- 
giant X-ray binaries (SGXBs), composed of a compact object 
orbiting around an early-type supergiant and fed by accretion 
from the strong radiative wind of the companion. These objects 
are persistent sources (Lx ~ 10 36 erg s _1 ), and their relative 
low number compared to the population of Be/X-ray binaries 
was explained as the consequence of the short lifetime of su- 
pergiant stars. 

The launch of the INTErnation al Gammay-Ray 



Astrophysics Laboratory (INTEGRAL, IWinkler et aU 12003) 
in October 2002 completely changed the situation, as many 
more HMXBs whose companion stars are supergiants were 
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discovered during the monitoring of the Galactic centre and 
the Galactic pl a ne usi n g the onboard IBIS/ISGRI instruments 
dUbertini et all 120031: Lebrunetai] 12003b . Most of these 



source s are reported in iBird et alj d2007l) and Bodagh ee et al 



(2007), and their studies have revealed two main features that 
were not present on previously known SGXBs: 

- first, many of them exhibit a considerable intrinsic absorp- 
tion, with a column density up to A^h ~ 2 x 10 24 cm~ 2 
in the case of IGR J 163 18-4848 dMatt & Guainazzill2003l) . 
which explains why previous high-energy missions had not 
detected them. 

- second, some of these new sources reveal a transitory na- 
ture. They are undetectable most of the time and occasion- 
ally present a fast X-Ray transient activity lasting a few 
hours. Moreover, they exhibit a quiescent luminosity of 
Lx ~ 10 33 erg s , well below the persistent state of other 
SGXBs. 



It then appears that the supergiant HMXBs discovered by 
INTEGRAL can be classified in two classes: one class of con- 
siderably obscured persistent sources that we will simply call 
obscured SGXBs in this paper and another of s upergiant fast 
X-ray transients (SFXTs. lNegueruela et"al]|2006bl) . 

High-energy observations can give some information about 
the compact object or about the processes that lead to the emis- 
sion but do not allow study of the companion star. It is there- 
fore very important to perform multi-wavelength observations 
of these sources - from optical-to-MIR wavelength - as this rep- 
resents the only way to characterise the companion or to detect 
dust around these highly obscured systems. However, positions 
given by INTEGRAL are not accurate enough (~ 2') to identify 
their optical counterparts, because of the large number of ob- 
jects in the error circle. Observations with X-ray telescopes like 
XMM-Newton or Chandra are therefore crucial because they 
allow a localisation with a position accuracy of 4" or better, 
which lowers the number of possible optical counterparts. 

We performed optical-to-MIR wavelength observations of 
several candidate SGXBs recently discovered with INTEGRAL. 
Optical and NIR observations were carried out at ESO/NTT us- 
ing EMMI and SofI instruments and aimed at constraining the 
spectral type of the companions through accurate astrometry, as 
well as the spectroscopy and photometry of the can didate coun 



terpar ts. They are reported in the companion paper ( Ch atv et al 



2008, CHA08 hereafter), and it is shown that most of these 



sources are actually supergiant stars. 

In this paper, we report MIR photometric observations of 
the companions of twelve INTEGRAL candidate SFXTs and 
obscured SGXBs that aimed at studying the circumstellar en- 
vironment of these highly absorbed sources and, more particu- 
larly, at detecting any MIR excess in their emission that could 
be due to the absorbing material. These sources were cho- 
sen because they had very accurate positions and confirmed 
2MASS counterparts. 

All the sources in our sample are listed in Table 1. The 
total galactic column density of neutral hydrogen A^h(Hi) is 
computed using the web version of th e A^h FTOOL from 
HEASARC. This tool uses the data from Dickey & Lockman 



( 1990), who performed H i observations from the Lyman-a' and 
21 cm lines. Moreover, A^h(H i) is the total galactic column den- 
sity, which means it is integrated along the line of sight over the 
whole Galaxy. Therefore, it is likely to be overestimated com- 
pared to the real value at the distance of the sources. 

The total galactic column density of molecular hydrogen 
AWH2) is computed us ing the velocity-integrated map (Wco ) 
and the X-ratio given in D ame. Hartmann. & Thaddeusld200lh . 
It is also likely to be overestimated compared to the real value 
at the distance of the sources, because it is integrated along 
the line of sight over the whole Galaxy. In contrast, A^hx, the 
intrinsic X-ray column density of the source, is computed 
from the fitting of the high-energy spectral energy distribution 
(SED), so it takes all the absorption into account at the right 
distance of the source. 

Using these observations, the results reported in the com- 
panion paper (CHA08), as well as archival photometric data 
from the USNO, 2MASS, and GLIMPSE catalogues when 
needed, we performed the broad-band SEDs of these sources 
and fitted them with a two-component black body model to as- 
sess the contribution of the star and the enshrouding material 
in the emission. The ESO observations, as well as our model, 
are described in Section 2. In Section 3, results of the fits for 
each source are given and these results are discussed in Section 
4. We conclude in Section 5. 



2. Observations 

2.1. MIR observations and data reduction 

The MIR observations were carried out on 2005 June 2 0-22 
and 2006 June 29-30 using VISIR dLagage etalJl2004 . the 
ESO/VLT mid-infrared imager and spectrograph, composed of 
an imager and a long-slit spectrometer covering several fil- 
ters in N and Q bands and mounted on Unit 3 of the VLT 
(Melipal). The standard "chopping and nodding" MIR observa- 
tional technique was used to suppress the background dominat- 
ing at these wavelengths. Secondary mirror-chopping was per- 
formed in the north-south direction with an amplitude of 16" 
at a frequency of 0.25 Hz. Nodding technique, needed to com- 
pensate for chopping residuals, was chosen as parallel to the 
chopping and applied using telescope offsets of 16". Because 
of the high thermal MIR background for ground-based obser- 
vations, the detector integration time was set to 16 ms. 

We performed broad-band photometry in 3 filters, PAH1 
(,i=8.59±0.42 urn), PAH2 (^=11.25+0.59 /mi), and Q2 
(/1=18. 72+0.88 //m) using the small field in all bands 
(1972xl9'.'2 and CC.'075 plate scale). All the observations were 
bracketed with standard star observations for flux calibration 
and PSF determination. The weather conditions were good and 
stable during the observations. 

Raw data were reduced using the IDL reduction package 
written by Eric Pantin. The elementary images were co-added 
in real-time to obtain chopping-corrected data, then the dif- 
ferent nodding positions were combined to form the final im- 
age. The VISIR detector is affected by stripes randomly trig- 
gered by some abnormal high-gain pixels. A dedicated destrip- 
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ing method was developed (Pantin 2008, in prep.) to suppress 
them. The MIR fluxes of all observed sources including the lcr 
errors are listed in Table 2. 



2.2. Archival data 

When we did not have optical-to-MIR data for our sources, we 
searched for the counterparts in 3 catalogues: 

- in the United States Naval Observatory (USNO) catalogues 
in B, R, and / for USNO-B1.0, B and R for USNO-A.2. 
Positions and fluxes accuracies are 0725 and 0.3 magni- 
tudes in the case of USNO-B.l, 072 and 0.5 magnitudes in 
the case of USNO-A.2. 

- in the 2 Micron All Sky Survey (2MASS), in J (1.25+0.16 
fim), H (1.65+0.25 pm) and Ks (2.17+0.26 pm) bands. 
Position accuracy is about 072. 

- in the Spitzer's Galactic Leg acy Infrared Mid-Plan e Survey 
Extraordinaire (GLIMPSE, Benjamin et al.l 120031) . survey 
of the Galactic plane Qb\ < 1° and between /=10° and 
1-65° on both sides of the Galactic centre) performed with 
the Spitzer Space Telescope, using the IRAC camera in four 
bands, 3.6+0.745 pm, 4.5+1.023 pm, 5.8+1.450 pm, and 
8+2.857 fim. 

All sources had a confirmed 2MASS counterpart and three of 
them (IGR J16195-4945, IGR J16207-5129, and IGR J16318- 
4848) had a GLIMPSE counterpart given in the literature. We 
found all the other GLIMPSE counterparts using the 2MASS 
positions and they are listed in Table 3. We used all the 
fluxes given in the GLIMPSE catalogue except in the case of 
IGR J 17252-36 16, IGR J 1739 1-3021 and IGR J 17544-26 19 
because their fluxes were not present in the catalogue tables. 
Nevertheless, we measured their fluxes on the archival images 
directly with aperture photometry. Uncertainties on the mea- 
surements were computed in the same way on the error maps 
given with the data. 



2.3. Absorption 

Absorption at wavelength A, A,j, is a crucial parameter to fit the 
SEDs, especially in the MIR. Indeed, inappropriate values can 
lead to a false assessment of the MIR excess. Visible absorption 
Ay was a free parameter of the fits. An accurate interstellar 

A 

absorption law - i.e. the wavelength dependence of the ratio 

Ay 

in the line of sight - was then needed to properly fit the SEDs. 

In the optical bands, we built the fu n ction with the ana- 
lytical expression given in Cardel li et al. (1989) who derived 



the average interstellar extinction law in the direction of the 

Galactic centre. Fro m 1 .25 pm to 8 yum, we used the analytical 

expression given in llndebetouw et aT ( 2005 ). They derived it 

from the measurements of the mean values of the colour excess 
(A -A i 

ratios .-' . K> from the colour distributions of observed stars 

(Aj-A K ) 

in the direction of the Galactic centre. They used archival data 
from 2MASS and GLIMPSE catalogues, which is relevant in 
our case as we use GLIMPSE fluxes. 



Above 8 pm, where absorption is dominated by the sil- 
icate features at 9.7 pm and 1 8 pm, we found several ex 
tincti o n laws in the litera ture dRieke et al. 1989; Lutz et al 



1996; Moneti et al. 2001), which exhibit some differences. 
Considering the high importance of a good assessment to cor- 

A 

rectly fit the MIR excess, we decided to assess the ratio -r- 1 - in 

Ay 

2 VISIR bands - PAH1 and PAH2 - from our data in order to 
bui ld the relevant law for our observations. 

iRieke & Lebofekvl dl985h gave the interstellar extinction 
law up to 13 pm, and from their results, we derived 0.043 < 

< 0.06. To get the best 



< 0.074 and 0.047 < 



Ay - ' - Ay 

values corresponding to our data in PAH1 and PAH2, we pro- 
ceeded in 3 steps. 

- We first selected the sources for which we had VISIR fluxes 
in PAH1 and/o r PAH2 and fit t ed th eir SEDs with extinction 



laws g iven in ICardelli et al.l (119891) and llndebetouw et al 



(2005) from 0.36 to 8 pm and half-interval values taken in 
PAH1 and PAH2. 
- Then, when we did not need any MIR excess to fit the IRAC 



fluxes, we adjusted the 



and 



ratios to improve 



Ay Ay 

the^f 2 of our fits. 
- We finally averaged all the extinction values obtained for 
all sources to get what we consider as the right ratios in 
PAH1 and PAH2 in the direction of the Galactic plane. 

The resulting values ar e in good agreeme nt with those given by 



the extinction law from lLutz et al. l[l996), so we chose their ex 
tinction law to fit our SEDs above 8 pm, the Q2 filter included. 

A 

The values we used in each band are listed in Table 4, and 

Ay 

the overall extinction law is displayed in Fig 1 . 



2.4. SEDs 

With all the archival and observational data from optical-to- 
MIR wavelength, we built the SEDs for these sources. We fit- 
ted them (using a x 2 minimisation) with a model combining 
two absorbed black bodies, one representing the companion 
star emission and a spherical one representing a possible MIR 
excess due to the absorbing material enshrouding the compan- 
ion star: 



AF(A) 



Inhc 1 
D t 2 A 4 



10 



-0.4A, 



R 2 



R» 2 



he 

,AkT, _ 1 



he 

,AkT u _ 1 



in W m 



-2 



We added to the flux uncertainties systematic errors as fol- 
lows: 

- a 2% systematic error in each IRAC band as given in the 
IRAC manuafl 

- comparing the variations of the flux calibration values ob- 
tained from standards with VISIR during our observation 
nights, we figured out that systematic errors with VISIR 
were about 5% at 10 pm and 10% at 20 pm. 



1 http://ssc.spitzer.caltech.edu/documents/som/ som8. 0.irac.pdf 
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The free parameters of the fits were the absorption in the 
V-band A v , the companion star black body temperature T„ and 
radius to distance ratio jy-, as well as the additional spherical 
component black body temperature and radius To and Rd- 



8 /J.m does not allow to reach definitive conclusions. Moreover, 
in both cases (with and without dust), the stellar component 
is consistent with a n O/B supergiant, as already suggested in 
Tomsick et ai1d2006l) . 



The best-fitting parameters for individual sources, as well as 2.2. IGR J1 6207-51 29 



corresponding x ®ce listed in Table 5 and the fitted SEDs 
are displayed in Fig 3. Moreover, 90%-confidence ranges of 
parameters are listed in Table 6. 

In Table 5, along with the best-fitting parameters, we also 
give the total galactic extinctions in magnitudes Am and Ah 2 in 
the line of sight, as well as the X-ray extinction of the source in 
magnitudes Ax- The values of Am, Am, and Ax are computed 
from A^h(Hi), Nnfjli), and jV » x given in Table 1 using the re 
lation Ah = 



3.1 



Rieke & Lebof skv 1985). 



-Nu feohlin. Savage. & Drakelll978 



3. Results 

3.1. IGR J 161 95-4945 

IGR J 16195-4945 was detected by INTEGRAL during obser- 
vatio ns carried out be tween 2003 February 27 and October 
19 JWalteretal.ll2004 and corresponds to the ASCA source 
AX J161929-4945 dSugizaki etafl |200lt ISidoli et alj Eool . 
Analysing INTEGRAL public data, Isidoli et all d20Q5l) derive 
an average flux level of ~ 17 mCrab (20-40 keV) . Performing 
a follow-up with INTEGRAL, ISguera et all d2006l) show it be- 
haves like an SFXT and report a peak-flux of ~ 35 mCrab (20- 

40 keV). 

Tom sick et al. d2006h observed the source with Chandra 



between 2005 April and July and give its position with 076 
accuracy. They fitted its high-energy emission with an ab- 
sorbed power law and derive F ~ 0.5 and A^h ~ 7 x 10 22 cirT 2 . 
Moreover, using their accurate localisation, they found its NIR 
and MIR counterparts in the 2MASS (2MASS J16193220- 
4944305) and in the GLIMPSE (G333.5571+00.3390) 
catalogues and performed its NIR photometry using ESO/NTT 
observations. They show its spectral type is compatible with 
an O, B, or A supergiant star. They also found possible 
USNO-A.2 and USNO-B.l counterparts. Nevertheless, as 
already suggested i n their paper, the USNO source is a blended 
foreground object dTovmassian et al.l l2006). 

We observed IGR J 16 195-4945 on 2006 June 30 in 
PAH1 during 1200 s, but did not detect it. Typical see- 
ing and airmass were 0788 and 1.07. We nevertheless fit- 
ted its SEP using the NIR and the GLIMPSE flux values 
given in Tom sick et al. ( 2006h and the best-fitting parameters 

R. 



are A v = 15.5, T,=23800 K, ^=5.96 x 10" 11 , 



T D =1160K, 

Rd=5.1R„, and the reduced x 2 is 3.9/2. 
The best-fitting parameters without the additional component 

are A v =16.1, T*=13800 K, ^=8.94 x 10", and the corre- 
sponding reduced^ 2 is 15.8/4. We then found a MIR excess. 

The additional component is then needed to correctly fit 
the SED, as this source exhibits a MIR excess. Nevertheless, as 
shown in Fig 4, this excess is small, and the lack of data above 



IGR J16207-5129 is an obscured SGXB that was discovered 
by INTEGRAL during obser vations carried out between 2003 
Fe bruary 27 and October 19 dWalter et al.ll2004 . 

Tomsick et al. d2006h observed it with Chandra during 
the same run as IGR J16 195-4945 and give its position 
with 076 accuracy. They also fitted its high-energy emis- 
sion with an absorbed power law and derive F ~ 0.5 and 
A^h ~ 3.7 x 10 22 cirT 2 . Thanks to their accurate locali- 
sation, they found its NIR and MIR counterparts in the 
2MASS (2MASS J 16204627-5 130060) and in the GLIMPSE 
(G333.4590+01.0501) catalogues and performed its NIR 
photometry using ESO/NTT observations. They show its 
temperature to be > 18000 K, which indicates the system is an 
HMXB. They also found its USNO-B1.0 (USNO-B1.0 0384- 
0560875) counterpart. 

We observed IGR J 16207-5 129 on 2006 June 29 in PAH1 
and PAH2 during 1200 s in each filter, and in Q2 during 2400 s. 
Typical seeing and airmass were 0772 and 1.09. We did not 
detect it in Q2 but in PAH1 and PAH2. The fluxes we derived 
are 21.7+1.4 mJy and 9.4+1.3 mJy, respectively. 

Using those values, as well as the fluxes from 
ESO/NT T observat ions and the GLIMPSE archives found in 



Tomsick et al 



2006), we fitted its SED, and the best-fitting 
parameters are A v =10.5, T» =33800 K, ]y=9.42 X 10~ n , 
and the reduced x 1 is 28.5/9. iNegueruela & Schurchl d2007l) 



find the spectral type is earlier than B1I; our parameters are 
therefore in good agreement with their results. 

The best fit with the additional component gives a larger 
reduced x 2 of 30/7 and To < 200 K, which is not significant, 
as the presence of such cold material marginally enhances the 
MIR flux. We therefore think IGR J 1 6207-5 1 29 is an O/B mas- 
sive star whose enshrouding material marginally contributes to 
its MIR emission. 



3.3. IGR J1 631 8-4848 

Main hi gh-energy characteristics of this source can b e 



found in Matt & Guainazzi (200J and IWalter et all 42003). 



IGR J16318-4848 was discovered by INTEGRAL on 2003 



January 29 dCourvoisier et alj I2003I) and was then observed 
with XMM-Newton, which allowed a 4" localisation. Those ob- 
servations showed that the source was exhibiting a strong ab- 
sorption of A^h ~ 2 x 10 24 cirT 2 , a temperature kT = 9 keV, and 
a photon index ~ 2. 



Using this accurate position, Filliatre & Chaty (2004) dis- 
covered its optical counterpart and co nfirmed the NIR coun- 
terpart proposed by IWalter et a"D d2003l) (2MASS J16314831- 
4849005). They also performed photometry and spectroscopy 
in optical and NIR on 2003 February 23-25 at ESO/NTT and 
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Table 1. Sample of sources studied in this paper. We give their name, their coordinates (J2000 and galactic), the total galactic 
column density of neutral hydrogen (/Vh(H I)) and the total galactic column density of molecular hydrogen (A^CH^)) in the line of 
sight, the X-ray column density of the source (/Vh x X their type (SFXT or OBS - obscured sources) and their spectral type (SpT). 
Their spe ctral classification s come fr om optical/NIR spec troscopy, reported in the follow ing r eferences (Ref): c : IChatv et alj 
d2008h. f: IFilliatre & Chatvl d2004. i: lin't Zand et all d2006l). nl: llSfegueruela et all d2005l) . n2: iNegueruela etal] d2006al) . n3: 
Nespoli et al. d2007h . p: iPellizza et al.l d2006h . t: iTomsick et al. I d2006h . 



Sources 


a(J2000) 


5(J2000) 


I 


b 


AWHDUO 22 


) A' H (H 2 )(10 22 


) A'nxdO 22 ) 


Type 


SpT 


Ref 


IGR J16195-4945 


16 19 32.20 


-49 44 30.7 


333.56 


0.339 


2.2 


5.4 


7.0 


SFXT? 


O/B 


t 


IGR J 16207-5 129 


16 20 46.26 


-51 30 06.0 


332.46 


-1.050 


1.7 


2.2 


3.7 


OBS 


O/B 


t 


IGR J 163 18-4848 


16 31 48.60 


-48 49 00.0 


335.62 


-0.448 


2.1 


3.6 


200.0 


OBS 


sgB[e] 


f 


IGR J 16320-4751 


16 32 01.90 


-47 52 27.0 


336.30 


0.169 


2.1 


4.4 


21.0 


OBS 


O/BI 


c 


IGR J 1635 8-4726 


16 35 53.80 


-47 25 41.1 


337.01 


-0.007 


2.2 


7.3 


33.0 


OBS 


sgB[e] 


c 


IGR J16418-4532 


16 41 51.00 


-45 32 25.0 


339.19 


0.489 


1.9 


3.6 


10.0 


SFXT? 


O/B 


c 


IGR J16465-4507 


16 46 35.50 


-45 07 04.0 


340.05 


0.135 


2.1 


5.9 


60.0 


SFXT 


B0.5I 


nl 


IGR J16479-4514 


16 48 06.60 


-45 12 08.0 


340.16 


-0.124 


2.1 


8.2 


7.7 


SFXT? 


O/BI 


c 


IGR J17252-3616 


17 25 11.40 


-36 16 58.6 


351.50 


-0.354 


1.6 


3.9 


15.0 


OBS 


O/BI 


c 


IGR J 1739 1-3021 


17 39 11.58 


-30 20 37.6 


358.07 


0.445 


1.4 


4.5 


30.0 


SFXT 


08Iab(f) 


n2 


IGR J 17544-26 19 


17 54 25.28 


-26 19 52.6 


3.26 


-0.336 


1.4 


7.8 


1.4 


SFXT 


091b 


P 


IGRJ19140+0951 


19 14 04.23 


+09 52 58.3 


44.30 


-0.469 


1.7 


3.7 


6.0 


OBS 


B1I 


n3 



Table 2. Summary of VISIR observations of newly discovered INTEGRAL sources. We give their MIR fluxes (mJy) in the PAH1 
(8.59 /mi), PAH2 (1 1.25 fim) and Q2 (18.72 fim) filters. When we did not detect a source, we give the upper limit. When no flux 
nor upper limit is given, the source was not observed in the considered filter. 



Sources 


PAH1 


PAH2 


Q2 


IGR J16195-4945 


< 6 


1 


< 7 


8 


< 50.3 


IGR J16207-5129 


21.7 ± 


1.4 


9.4 ± 


1.3 


< 53.4 


IGR J16318-4848 


426.2 ± 3.0 


317.4 + 3.4 


180.7+15.3 


IGR J 16320-4751 


12.1 + 


1.7 


6.3 ± 


1.8 




IGR J16358-4726 


< 6 


9 








IGR J16418-4532 


< 5 


8 








IGR J 16465-4507 


6.9 ± 


1.1 


< 5 







IGR J 16479-45 14 


10.9 + 


1.2 


7.0 ± 


1.6 




IGR J17252-3616 


6.1 ± 


0.6 


< 5 







IGR J 1739 1-3021 


70.2 + 


1.6 


46.5 + 


2.6 




IGR J 17544-26 19 


46.1 + 


2.8 


20.2 + 


2.1 




IGRJ19140+0951 


35.2 + 


1.4 


19.1 + 


1.4 





show that the source presents a significant NIR excess and that 
it is strongly absorbed (Ay ~ 17.4). The spectroscopy revealed 
an unusual spectrum with a continuum very rich in strong emis- 
sion lines, which, together with the presence of fo rbidden lines, 
points towards an sgB fel companion star (see e.g. lLamers et al. 
1998: Zickgradl999l for definition and characteristics of these 
stars). 

Using the 2MASS magnitud es, the GLIMPSE (G3 35.6260- 
00.4477), and the MSX fluxes, iKaplan et alJ d2006h fitted its 
SED with a combination of a stellar and a dust component 
black bodies, and shows that the presence of warm dust around 
the system was necessary for explaining the NIR and MIR 
excess. From their fit, they derive Ay ~ 18.5, Td=1030 K, and 
R D =10R t . 



We observed IGR J16318-4848 with VISIR twice: 



- the first time on 2005 June 21 during 300 s in PAH1 and 
PAH2, and 600 s in Q2. Typical seeing and airmass were 
0781 and 1.14. We detected the source in all bands, and 
the derived fluxes are 409.2+2.4 mJy, 322.4+3.3 mJy, and 
172.1 + 14.9 mJy in PAH1, PAH2, and Q2, respectively. 

- the second on 2006 June 30 during 600 s in all bands. 
Typical seeing and airmass were 0768 and 1.09. We de- 
tected the source in all bands, and the derived fluxes are 
426.2+3.0 mJy, 317.4+3.4 mJy, and 180.7+15.3 mJy in 
PAH1, PAH2, and Q2, respectively. 

Those observations show that IGR J16318-4848 is very bright 
in the MIR (it is actually the brightest source in our sample) 
and that its flux was constant within a year, considering VISIR 
systematic errors. 

Usin g data from our last run , as well as the magnitudes 



given in IFilliatre & Chatvl (|2004) in the optical and the NIR, 



and the fluxes from the GLIMPSE archives, we fitted its 
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Table 3. List of GLIMPSE counterparts we found for 9 sources. We give their name, their separation from the 2MASS counter- 
parts and their fluxes in mJy. 



Sources 


GLIMPSE counterpart Separation 


3.6yU 


m 


4.5 yum 


5.8/jm 


8/jm 


IGR J 16320-4751 


G336.3293+00.1689 


«' 17 


48.2 + 


1.9 


44.3 + 2.1 


36.0 ± 2.0 


17.3 ± 1.0 


IGR J16358-4726 


G337.0994-00.0062 


0T46 


5.9 + 0.5 


5.6 + 0.6 


5.3 ± 1.7 




IGR J16418-4532 


G339.1889+00.4889 


«' 28 


12.5 ± 


0.9 


9.5 + 0.6 


5.6 + 0.6 


3.6 ± 0.4 


IGR J 16465-4507 


G340.0536+00.1350 


a/ 16 


45.0 ± 


2.0 


32.6+1.5 


22.0 ± 0.9 


13.7 + 0.6 


IGR J 16479-45 14 


G340. 1630-00. 1239 


a' 13 


68.6 ± 


3.1 


49.6 ± 2.0 


41.2 + 2.3 


19.4 + 0.9 


IGR J 17252-36 16 






32.6 ± 


3.7 


24.7 ± 4.7 


21.8 + 5.0 


9.6 + 6.5 


IGR J 1739 1-3021 






375.9 ± 


44.0 


297.0 + 31.0 


205.0 ± 33.0 


111.0 + 28.0 


IGR J 17544-26 19 






213.9 + 


25.4 


137.0+18.9 


99.6 + 7.6 


66.5 ± 12.1 


IGR J19140+0951 


G044.2963-00.4688 


0748 


185.0 4 


:9.3 


152.0+11.1 


103.9 + 5.4 


62.0 + 2.1 



Table 4. Adopted ^ values. 



Filters 


U 


B V R I J H Ks 3.6yum 4.5yum 5.8yum 


8 yum 


8.59 (im 


11.25 /mi 


18.72 (im 


A, 

A,. 


1.575 


1.332 1 0.757 0.486 0.289 0.174 0.115 0.0638 0.0539 0.0474 


0.0444 


0.0595 


0.0605 


0.040 



10 



< 

< 



0.1 



0.01 



CardeNi et al. 


Indebetouw et al. 


Lutz et al. 









0.1 



10 



100 



X ((j.m) 



Fig. 1. Adopted extinction law. We used the la w given in ICardelli et al. ( 1989b in the optical, the one given in llndebetouw et al 
d2005l) from 1.25 fim to 8 /im, and the law from lLutz et al.1 (1 19961) above 8 fan. 



SED and the best-fitting parameters are Ay = 17, T» =22200 K, 
^=3.74xlO" 10 , T D =1100 K, R d =10R», and the reduced^ 2 is 
6.6/6. The best-fitting parameters without the additional com- 
ponent are A v =17.9, T»=18200 K, ]y-=5.1 x 10" 10 , and the 

corresponding reduced^- 2 is 425/8. We then confirm that the 
MIR excess is likely due to the prese nce of warm dust around 
the system, as already suggested b y iFilliatre & Chatvl (|2004) 
and reported in Kapl an et alj d2006h . 



3.4. IGR J1 6320-4751 

IGR J 16 320-4751 was detec ted by INTEGRAL on 2003 
February dTomsick et al. 2003b and corresponds to the ASCA 
source AX J 163 1.9-4752. iRodriguez et al] d2003b report ob- 
servations with XMM-Newton. They give an accurate localisa- 



tion (3") and fitted its high-energy spectrum with an absorbed 
po wer law. They derive T ~ 1.6 andA^H ~ 2.1 x 10 23 cm" 2 . 



Lutovinov et al.l (12005b report the discovery of X-Ray pul- 
sations (P~ 1309 s), which proves the c ompact object is a 



neutron star. Moreover, ICorbet et al.l {2004) obtained the light 
curve of IGR J 16320-4751 between 2004 December 21 and 
2005 September 17 with Swift and report the discovery of a 
8.96 days orbital period. IGR J16320-4751 is then an X-ray 
bin ary whose compact object is a neutron star. 



Neg ueruela & Schurchl (I2007h searched for the NIR coun 



terpart of the source in the 2MASS catalogue and found its po- 
sition was consistent with 2MASS J16320215-4752289. They 
also concluded that, if it was an O/B supergiant, it had to be 
extremely absorbed. 

The optical and NIR photometry and spectroscopy of this 
source were carried out at ESO/NTT, and results are reported 
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Table 5. Summary of best-fitting parameters of the SEDs of the sources. We give the total galactic extinctions in magnitudes A Hl 
and A H ,, the X-ray extinction of the source in magnitudes A x and then the parameters themselves: the extinction in the optical 

Ay, the temperature T* and the jj- ratio of the companion and the temperature and radius Td and Rd (in R* unit) of the dust 

component when needed. We also add the reduced^ 2 we reach for each fit. 



Sources 


A H i 


A H2 


A x 


A v 


T.(K) 


R. 
D, 




T D (K) R 


d(R.) 


* 2 /dof 


IGR J16195-4945 


11.7 


28.6 


37.4 


15.5 


23800 


5.96 x 10" 


ii 


1160 


5.1 


3.9/2 


IGR J 16207-5 129 


9.3 


11.8 


19.8 


10.5 


33800 


9.42 x 10" 


n 






28.5/9 


IGR J 163 18-4848 


11.0 


19.2 


1069.5 


17.0 


22200 


3.74 x 10" 


10 


1100 


10.0 


6.6/6 


IGR J 16320-4751 


11.4 


23.6 


112.2 


35.4 


33000 


1.38 x 10" 


in 






7.7/6 


IGR J 1635 8-4726 


11.8 


39.4 


176.4 


17.6 


24500 


3.16 x 10- 


1 1 


810 


10.1 


3.6/2 


IGR J16418-4532 


10.1 


19.3 


53.5 


14.5 


32800 


3.77 x 10" 


n 






1.4/4 


IGR J16465-4507 


11.3 


31.4 


320.7 


5.9 


25000 


6.40 x 10- 


1 1 






13.9/7 


IGR J16479-4514 


11.4 


43.8 


41.2 


18.5 


32800 


1.00 x 10- 


in 






7.4/6 


IGR J17252-3616 


8.3 


20.1 


80.2 


20.8 


32600 


7.57 x 10" 


n 






3.8/5 


IGR J 1739 1-3021 


7.3 


23.9 


160.4 


9.2 


31400 


1.80 x 10" 


in 






11.7/10 


IGR J 17544-26 19 


7.7 


41.5 


7.70 


6.1 


31000 


1.27 x 10- 


in 






6.1/8 


IGRJ19140+0951 


9.0 


14.0 


32.1 


16.5 


22500 


1.92 x IO- 


L0 






14.4/6 



Table 6. Ranges of parameters that give acceptable fits (90%-confidence) for each source. 



Sources 


AA 


V 


AT. 




A 


R » 
D, 




AT D 


ARd 


IGR J16195-4945 


14.8- 


15.8 


13100 


- 25900 


5.68 x 10" 


1 1 


- 6.68 x 10" 


ii 


950- 


1460 


3.9 - 6.6 


IGR J 16207-5 129 


10.4- 


10.6 


25200 


- 36000 


9.36 x 10- 


n 


- 1.05 x 10" 


10 








IGR J 163 18-4848 


16.6- 


17.4 


19300 


- 24500 


3.65 x 10" 


10 


- 3.84 x 10" 


10 


960- 


1260 


8.8- 11.8 


IGR J 16320-4751 


34.8- 


35.5 


22000 


- 35600 


1.33 x 10- 


10 


- 1.69 x 10" 


10 








IGR J 1635 8-4726 


17.1 - 


18.1 


19500 


- 36000 


2.64 x 10" 


11 


-3.52x 10" 


11 


630- 


1020 


8.0-13.8 


IGR J16418-4532 


13.6- 


14.7 


10600 


- 36000 


3.58 x 10- 


1 1 


- 5.44 x 10" 


1 1 








IGR J16465-4507 


5.0- 


6.1 


15400 


- 33600 


5.50 x 10- 


1 1 


- 7.95 x 10" 


1 1 








IGR J16479-4514 


18.4- 


18.8 


26200 


- 36000 


9.48 x 10" 


II 


- 1.13 x 10" 


10 








IGR J17252-3616 


20.3- 


21.0 


20500 


- 36000 


7.17 x 10" 


11 


- 9.48 x 10" 


1 1 








IGR J 1739 1-3021 


8.8- 


9.4 


16100 


- 32200 


1.78 x 10" 


10 


- 2.63 x 10" 


10 








IGR J 17544-26 19 


6.0- 


6.2 


26700 


- 35500 


1.18 x 10" 


10 


- 1.38 x 10" 


10 








IGRJ19140+0951 


15.7- 


16.7 


13200 


-28100 


1.71 x 10" 


10 


-2.53x IO- 


10 









in CHA08. It is shown that its NIR spectrum is consistent 
with an O/B supergiant and that its intrinsic absorption is very 
high, because it was not detected in any of the visible bands. 
We searched for the MIR counterpart of IGR J 16320-4751 
in the GLIMPSE archives and found it to be consistent with 
G336.3293+00.1689. 

We observed IGR J 16320-4751 with VISIR on 2005 June 
20 in PAH1 and PAH2, and the respective exposure times were 
1800 s and 2400 s. Typical seeing and airmass were Of.' 63 and 
1.13. We detected it in both filters, and the respective fluxes 
are 12.1 + 1.7 mJy and 6.3+1.8 mJy. Using the ESO/NTT NIR 
magnitudes given in CHA08, as well as the GLIMPSE and the 
VISIR fluxes, we fitted its SED and the best-fitting parameters 
are A v =35.4, T„=33000 K, ^ = 1.38 x 10 10 , and the reduced 

X 2 is 7.7/6. This result is in good agreement with an extremely 
absorbed O/B supergiant as reported in CHA08. 

The best fit with the additional component gives a larger 
reduced x 2 of 8/4 and T D < 200 K. We therefore think that 



IGR J 1 6320-475 1 is an O/B supergiant whose enshrouding ma- 
terial marginally contributes to its MIR emission, even if its 
intrinsic absorption is extremely high. 

3.5. IGR J1 6358-4726 

IGR J1635 8-4726 was detected w ith INTEGRAL on 2003 
March 19 ( Revnivtsev et al. 2003b and first observed with 



Chandra on 2003 March 24 dPatel et all 120041) . They give its 



position with (X'6 accuracy and fitted its high-energy spec- 
trum with an absorbed power law. They derive F ~ 0.5 and 
Nh ~ 3.3 x 10 23 cm" 2 . They also found a 5880+50 s mod- 
ulation, which could be either a neutron star p ulsation or an 
orbital modulation. Nevertheless, Pate l et al. (2006) performed 
detailed spectral and timing analysis of this source using multi- 
satellite archival observations and identified a 94 s spin up, 
which points to a neutron star origin. Assuming that this spin 
up was due to accretion, they estimate the source magnetic field 
is between 10 13 and 10 15 G, which could support a magnetar 
nature for IGR J16358-4726. 
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Kouveliotou et all d2003l) propose 2MASS J16355369- 
4725398 as the possible NIR counterpart, and NIR spec- 
troscopy and photometry of this counterpart was performed at 
ESO/NTT and is reported in CHA08. They show that its spec- 
trum is consistent with a B supergiant belonging to the same 
family as IGR J16318-4848, the so-called B[e] supergiants. 
We also found its MIR counterpart in the GLIMPSE archives 
(G337.0994-00.0062). 

We observed IGR J16358-4726 with VISIR on 2006 June 
29 but did not detect it in any filter. Using the NIR magni- 
tudes given in CHA08 and the GLIMPSE fluxes, we fitted its 
SED and the best-fitting parameters are Ay=17.6, T„=24500 K, 

^ = 3 . 1 6 x 1 0" 1 1 , T D = 8 1 K, R D = 1 0. 1 R„ and the reduced X 2 is 
3.6/2. The best-fitting parameters without the additional com- 
ponent are A v = 16.7, T»=9800 K, ]^=6.05 x 10" and the 
corresponding reduced^ 2 is 8.8/4. 

The additional component is then necessary to correctly fit 
the SED, since this source exhibits a MIR excess (see Fig 4). 
Even if we lack MIR data above 5.8 fim, we think this excess 
is real and stems from warm dust, as it is consistent with the 
source being a sgB[e], as reported in CHA08. 



Observations were carried out on 2004 Septembe r 14 with 
XMM-Newton, and IZurita Heras & Walter! d2004 report a 
position with 4" accuracy, allowing identification of an NIR 
counterpart in the 2MASS catalogue (2MASS J16463526- 
4507045 =USNO - B 1.0 448-00520455). With the ESO/NTT, 
Negueruela et al.l (120051) performed intermediate-resolution 



3.6. IGR J 1641 8-4532 

IGR J16418-4532 was 
2003 Februar y 1-5 (ITomsick et al 

observations, Sguera et al.l (2006) report an SFXT be- 
haviour of this source and a peak-flux of ~ 80 mCrab 
(20-30 keV). Moreover, us i ng XM M-Newton and INTEGRAL 



discovered w ith INTEGRAL on 
20041) . Using INTEGRAL 



observations, Walter et al 



(2006) report a pulse period of 
1246+100 s and derive Nn ~ 10 23 crrT 2 . They also proposed 
2MASS J16415078-4532253 as its likely NIR counterpart. 
The NIR photometry of this counterpart was performed at 
ESO/NTT and is reported in CHA08. We also found the MIR 

counterpart in the GLIMPSE archives (G339.1889+004889). 3.5. IQR J1 6479-451 4 



spectroscopy of the source, estimate the spectral type is a 
B0.5I, and propo se that it is a n SFX T. Using XMM-Newton 
and INTEGRAL, IWalter et all d2006l) find a pulse period of 
227+5 s and derive iV H ~ 6 x 10 23 cm" 2 . We found its MIR 
counterpart in the GLIMPSE archives (G340.0536+00.1350) 
using the 2MASS position. 

We observed IGR J16465-4507 with VISIR twice: 

- The first one on 2005 June 20 during 600 s in PAH1. 
Typical seeing and airmass were 0781 and 1.14. We de- 
tected the source and the derived flux is 8.7+1.8 mJy. 

- the second on 2006 June 30 during 1200 s in PAH1 and 
PAH2. Typical seeing and airmass were 0768 and 1.09. We 
detected the source in PAH1 but not in PAH2. The derived 
flux is 6.9+1.1 mJy. 

These observations show that the IGR J16465-4507 MIR flux 
was constant during the year. 

Using the USNO-B1.0, 2MASS, and GLIMPSE flux val- 
ues, as well as our VISIR data, we fitted its SED, and the best- 
fitting parameters are A v =5.9, T»=25000 K, ]j=6.4 x 10 -11 , 

and the reduced^- 2 is 13.9/7. The best fit with the additional 
components gives a larger reduced x 2 of 20/5 and To < 200 
K. 

We then conclude that no additional component is needed 
to explain the MIR emission of this source, and the parameters 
derived from our fit are in good agreement with IGR J 16 465- 
4507 being a B0. 51 as reported in lNegueruela et al.l (|2005). 



We observed IGR J16418-4532 with VISIR on 2006 June 
29 but did not detect it in any filter. Using The NIR magni- 
tudes given in CHA08, as well as the GLIMPSE fluxes, we 
fitted its SED, and the best-fitting parameters are Av=14.5, 
T t =32800 K, ]^=3.77 x 10 -11 , and the reduced^ 2 is 1.4/4. 
The best fit with the additional component gives a larger re- 
duced x 1 of 3.9/2 and T D < 200 K. 

Uncertainties on the data are high, which is the reason why 
the reduced x 2 are low. Nevertheless, parameters of the fit, as 
well as the 90%-confidence ranges of parameters listed in Table 
6, are consistent with an O/B massive star nature. The temper- 
ature of the additional component being insignificant, we con- 
clude this source is an O/B massive star whose enshrouding 
material marginally contributes to its MIR emission. 

3.7. IGR J 16465-4507 

IGR J16465-4507 is a transient source discovered with 



IGR J16479-4514 was discovered with INTEGRAL on 2003 



August 8-9 dMolkov et al.ll2003l) . ISguera et al l (12005I) suggest 
it is a fast transient aft er they detected recurrent outbursts, 
and ISguera et aL ( 2006I) re p ort a peak-flux of ~ 120 mCrab 



INTEGRAL on 2004 September 6-7 (Lu tovinovet all 12004). 



(20-60 keV). IWalter et al.l (120061) observed it with XMM- 
Newton and gave its position with 4" accuracy. Moreover, 
they derive Nh ~ 7.7 x 10 22 crrT 2 from their observations. 
They also propose 2MASS J16480656-4512068=USNO- 
B1.0 0447-0531332 as its likely NIR counterpart. The NIR 
spectroscopy and photometry of this counterpart were per- 
formed at ESO/NTT and are reported in CHA08. It is shown 
that its spectrum is consistent with an O/B supergiant. We 
also found the MIR counterpart in the GLIMPSE archive 
(G339. 1889+004889). 

We observed IGR J16479-4514 with VISIR on 2006 June 
29 in PAH1 and PAH2, and the exposure time was 1200 s 
in each filter. Typical seeing and airmass were 079 and 1.14. 
We detected it in both filters, and the respective fluxes are 
10.9+1.2 mJy and 7.0+1.6 mJy. Using the NIR magnitudes 



F. Rahoui et al.: MIR observations revealing the obscured INTEGRAL binary systems 



9 



given in CHA08, as well as the GLIMPSE and the VISIR 
fluxes, we fitted its SED, and the best-fitting parameters are 

A v = 1 8 .5, T. =32800 K, ^ = 1 .00 x 1 0" 10 , and the reduced X 2 is 
7.4/6. The best fit with the additional component gives a larger 
reduced x 1 of 9/4 and T D < 200 K. 

We then do not need any additional component to fit the 
SED, and our result is consistent with IGR J16479-4514being 
an obscured O/B supergiant, in good agreement with CHA08. 

3.9. IGR J1 7252-361 6 

IGR J 17252-36 16 is a heavily-absorbed persistent source dis- 
covered with INTE GRAL on 2004 February 9 and reported in 



also give its optical/NIR counterpart 2MASS J17391 155- 
3020380=USNO-B1.0 0596-0585865 and classify IGR 



Walter et alJ (120041) . It was observed w i th XM M-Newton on 
2004 March 21, and IZurita Heras et al.1 (|2006) give its posi- 
tion with 4" accuracy. Using the XMM-Newton observations, as 
well as those carried out with INTEGRAL, they show the source 
was a binary X-ray pulsar with a spin period of ~ 413.7 s and an 
orbital period of ~ 9.72 days, and derive Nh ~ 1.5 X 10 23 crrT 2 . 
Moreover, they fitted its high-energy spectrum with either an 
absorbed compton (kT~ 5.5 keV and t ~ 7.8) or a flat power 
law (F ~ 0.02). 

In their paper, they propose 2MA SS J 17251 139-3616575 
as its likely NIR counterpart, as do Negueruela & Schurchl 
(2007). The NIR spectroscopy and photometry of this counter- 
part were performed at ESO/NTT and are reported in CHA08, 
where it is shown that its spectrum is consistent with an O/B 
supergiant. Using the 2MASS position, we searched for its 
MIR counterpart in the GLIMPSE catalogue. Unfortunately, 
we did not find its IRAC fluxes in the database. Nevertheless, 
we found post-Basic Calibrated Data (post-BCD) images of 
the source in all filters. We then reduced those data and de- 
rived fluxes directly from the images. They are listed in Table 3. 

We observed IGR J 17252-36 16 with VISIR on 2006 June 
30 in PAH1 and PAH2 and the exposure time was 1200 s in 
each filter. Typical seeing and airmass were 0797 and 1.09. We 
detected it in PAH1, and the derived flux is 6.1+0.6 mJy. Using 
the NIR magnitudes given in CHA08, as well as the GLIMPSE 
and the VISIR fluxes, we fitted its SED, and the best-fitting pa- 
rameters are A v =20.8, T. =32600 K, ^=7.57 x lO -11 , and the 

reduced* 2 is 3.8/5. The best fit with the additional component 
gives a larger reduced* 2 of 6.9/3 and Tp < 200 K. We then 
do not need any additional component to fit the SED, and our 
result is consistent with IGR J 17252-36 16 to be an obscured 
O/B supergiant, in good agreement with CHA08. 

3.10. IGR J1 7391 -3021 

IGR J17391-3021 is a transient source discovered with 
INTEGRAL on 2003 August 26 dSunvaev et alJ I2003I) and 



it corresponds to the Rossi X-ray Timing Explorer (RXTE) 
source XTE J 1 739-302. ISguera et al.1 ( 12005b analysed archival 
INTEGRAL data and classified the source as a fast X- 
ray transient pre senting a typical neutron star spectrum. 
Smith et al. I d2006l) observed it with Chandra on 2003 October 
15 and give its precise position with 1" accuracy. They 



J17391-3021 as an SFXT. lNegueruela et al.1 (l2006al) performed 
optical/NIR photometry and spectroscopy of the companion 
using ESO/NTT and find it is a 08Iab(f) star whose distance is 
~ 2.3 kpc. CHA08 also report optical and NIR spectroscopy 
and photometry of the companion carried out at ESO/NTT 
and confirm the nature of the companion. Using the 2MASS 
position, we searched for its MIR counterpart in the GLIMPSE 
catalogue and as for IGR J17252-3616, we had to reduce 
post-BCD data and derive the fluxes directly from the images. 
The fluxes are listed in Table 3. 

We observed IGR J17391-3021 with VISIR on 2005 June 
20 in PAH1 and PAH2, and the exposure time was 600 s in each 
filter. Typical seeing and airmass were 0763 and 1.13. We de- 
tected it in both filters, and the derived fluxes are 70.2+ 1 .6 mJy 
and 46.5+2.6 mJy. Using the optical and the NIR magnitudes 
given in CHA08, as well as the GLIMPSE and the VISIR 
fluxes, we fitted its SED and the best-fitting parameters are 
A v =9.2, T„=31400 K, ]j = 1.8 x 10~ 10 , and the reduced* 2 
is 11.7/10. The best fit with the additional component gives a 
larger reduced * 2 of 1 5 . 3/8 and T D < 200 K. 

We then do not need any additional component to fit the 
SED, and the parameters derived from our fit are in good agree- 
ment with IGR J 1739 1-3021 to be an Q8Iab(f) supergiant star, 
as initially reported in lNegueruela et al. I d2006ah . 



3.11. IGR J1 7544-261 9 

IGR J 17544-26 19 is a transient source discovered with 



INTEGRAL on 2003 Septem ber 17 dSunvaev et alJ 120031) . 
Gonzalez-Riestra et all d2004h observed it with XMM-Newton 
and derive Nh ~ 2 x 10 22 crrT 2 . They also confirm the 
association of the companion with 2MASS J17542527- 
2619526=USNO- B 1.0 0636 - 06209 33, as proposed in 
Rodriguez d2003l) . lin'tZandl d2005l) report on observations 
performed with Chandra, give its position with 076 accuracy, 
A'h ~ 1.36xl0 22 crrT 2 , and show that its high-energy spectrum 
is typical of an accreting neutron st ar. Moreover, they i dentify 
the counterpart as a blue supergiant. lSguera et al. d2006l) report 
a peak -flux of ~ 240 mCrab. Using ESO/NTT IPellizza et al ' 



(2006) performed optical/NIR spectroscopy and photometry 
of the companion, and give its spectral type as 091b at 
2.1-4.2 kpc. Using the 2MASS position, we searched for 
its MIR counterpart in the GLIMPSE catalogue, and as for 
IGR J17252-3616 and IGR J17391-3021, we had to reduce 
post-BCD data and derived the fluxes directly from the images. 
The fluxes are listed in Table 3. 

We observed IGR J 17544-26 19 with VISIR on 2005 June 
20 in PAH1 and PAH2, and the exposure time was 600 s in 
PAH1 and 1200 s in PAH2. Typical seeing and airmass were 
0764 and 1.13. We detected it in both filters, and the derived 
fluxes are 46.1+2.8 mJy and 20 .2+2.1 mJy. Using the magni- 
tudes from lPellizza et al.1 d2006|), the GLIMPSE and the VISIR 
fluxes, we fitted its SED, and the best-fitting parameters are 
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A v =6.1, T, =31000 K, ]^ = 1.27 x 10~ 10 , and the reduced* 2 is 
6.1/8. The best fit with the additional component gives a larger 
reduced * 2 of 9/6 and T D < 200 K. 

We then do not need any additional component to fit the 
SED, and the parameters derived from our fit are in good agree- 
ment with IGR J 17544-26 19 as an Q9 Ib supergiant star, as ini- 
tially reported in lPellizza et al. I d2006h . 




3.12. IGR J1 91 40+0951 



IGR J19140+0951 is a persistent source t hat was discovered 



with INTEGRAL on 2003 March 6-7 (lHannikainen et al 
2003D . Observations carried out with RXT E allowed T ~ 1.6 
and jV H ~ 6 x 10 22 cirT 2 to be derived dSwank & Markwardt 
20031) . Timing analysis of the R XTE data showed a period 
of 13.55 days (iCorbet et alJ 120041) . which shows the binary 
nature of the source. After a comprehensiv e analy sis of 
INTEGRAL and RXTE data, iRodriguez etal] d2005b show 
the source is spending most of its time in a faint state but 
report high variations in luminosity and absorption column 
density (up to ~ 10 23 cm -2 ). They also find evidence that 
the compact object is a neutron star rather than a black 
hole. Using Chandra observations carried out 2004 May 11, 
in't Zand et al. d2006l) give its position with 076 accuracy. 
This allowed them to find its NIR counterpart in the 2MASS 
catalogue (2MASS J19140422+0952577). Moreover, they 
searched for its M IR counterp art in the Mid-course Space 
Experiment (MSX iMilll 1 1994 and found an object at 8.3 
yum. The NIR photometry and spectroscopy of this source 
were performed at ESO/NTT and results are reported in 
CHA08. It is shown that its spectrum is consistent with an O/B 
massive star, in good agreement with iNespoli et alJ (|2007l) . 
who show it is a B1I supergiant. Using the 2MASS position, 
we also found its MIR counterpart in the GLIMPSE archive 
(G044.2963-00.4688). 

We observed IGR J19140+0951 with VISIR on 2006 June 
30 in PAH1 and PAH2, and the exposure time was 1200 s in 
each filter. Typical seeing and airmass were l'/12 and 1.17. 
We detected it in both filters, and the derived fluxes are 
35.2+1.4 mJy and 19.1 + 1.4 mJy. We point out that the ob- 
ject given as the MS X counterpart of IGR J19140+0951 in 
in t Zand et alJ (2006) is a blended source. Indeed, VISIR im- 
ages, whose resolution is far better, clearly show there are two 
sources in the field, IGR J19140+0951 and a very bright south- 
ern source (see Fig 2). 

Using the magnitudes given in CHA08, as well as the 
fluxes from GLIMPSE and our observations with VISIR, we 
fitted its SED, and the best-fitting parameters are Av=16.5, 
T t =22500 K, ]^ = 1.92 x 10 -10 , and the reduced X 2 is 14.4/6. 
The best fit with the additional component gives a larger re- 
duced* 2 of 20.2/4 and T D < 200 K. 

We do not need any additional component to fit the SED, 
and the parameters derived from our fit are in good agreement 
with IGR J 19 140+0951 to be an B II supergiant star, as initially 
reported in lNespoli et alJ ( 12007b . 




Fig. 2. VISIR image of IGR J19140+0951 in PAH1 (8.59 ^m). 
19'.'2xl9'.'2 field of view and 07075 plate scale. We clearly see 
the two sources that were blended with MSX. The MIR coun- 
terpart of IGR J19140+0951 is the northern source. 

4. Discussion 

4.1. Bvs B[e] supergiants stellar winds 

All SEDs were best-fitted without any dust component (even 
the very absorbed one like IGR J16320-4751), except three of 
them (IGR J16195-4945, IGR J163 18-4848, and IGR J16358- 
4726, see Fig. 3) that exhibit a MIR excess likely due to the 
presence of dust in their stellar wind. 

Blue supergiants are known to exhibit a very strong but 
sparse stellar wind of high velocity (~ 1000-2000 km s~'). 
This has been explained through the so-called radi ation 
line-driven CAK model dCastor. Abbott. & Kleinl 1 19751) in 
which the wind is driven by absorption in spectral lines. Hot 
stars emit most of their radiation in the ultraviolet (UV) where 
their atmosphere has many absorption lines. Photons coming 
from the photosphere of the star with the same wavelength are 
absorbed and re-emitted to the expanding medium in a random 
direction with almost the same momentum, which results 
in acceleration of the wind. This process is very effective 
because the line spectrum of the scattering ions in the wind 
is Doppler-shifted compared to the stellar rest frame, so the 
scattering atoms are shifted with respect to their neighbours at 
lower velocities and can interact with an unaffected part of the 
stellar spectrum. 

IGR J16318-4848 was proven to belong to a particular 
class of Bl supergian ts, the B[e] supergiants or sgB[e] 
dFilliatre&Chatyll2004. A physical definition of B[e] stars 
can be found in lLamers et al.l (11998). We just recall two of 
the characteristics here: the presence of forbidden emission 
lines of [Fen] and [Oi] in the NIR spectrum and of a strong 
MIR excess due to hot circumstellar dust that re-emits the 
absorbed stellar radiation through free-free emission. An 
sgB[e] is defined by the B[e] phenomenon, the indication of 
mass-loss in the optical spectrum (P-cygni profiles), and a 
hybrid spectrum characterised by the simultaneous presence 
of narrow low-excitation lines and broad absorption features 
of high-excitation lines. This hybrid nature was empirically 
explained by the simultaneous presence of a normal supergiant 
hot polar wind (fast and sparse) and responsible for the broad 
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Fig. 3. Optical-to-MIR absorbed (line) and unabsorbed (dotted-line) SEDs of 12 INTEGRAL sources, including broad-band photometric data 
from ESO/NTT, 2MASS, GLIMPSE, and VISIR. 



lines and a cool equatorial outflowing disk- like wind (slow and 
dense) responsible for the narrow line s ( Shore & Sanduleak 
19831; IZick?raf et alJll98H: IShore et alJll987b . This empirical 
model has received som e confirmation from polarimetry 
dOudmaijer & Drewlll999h . 



Kraus & M iroshnichenko (2006) 



There are a few models that explain the creation of this 
outflowing disk, and all of them consider the star rotation to 
be an important parameter in the process. In this paper, we 
present only the most consistent of them, the Rotation Induced 
Bi-stability mechanism (RIB), but a review can be found in 



The lines responsible for the creation of the wind are depen- 
dent on the ionisation structure, and a change in this structure 
leads to a cha nge in the radiative flux. Th is is the bi-stability 
jump found bv lLamers & Pauldrachldl99ll) . which appears for 
B stars with effective temperatures of about 23000 K. Above 
this temperature, the wind tends to be fast and sparse. Below, 
the mass-loss rate if five times higher and the terminal velocity 
two times slower, which leads to a wind that is ten times denser. 
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Fig. 4. IR SEDs of IGR J16195-4945, IGR J16318-4848, and 
IGR J16358-4726 in the NIR and the MIR. We show their 
SEDs including the contribution of the star and the dust (line), 
the star only (dashed-line), and the dust only (dotted-line). 



Cassinelli & Ignacd (I1997I) propose that the same effect is 



important from polar to equatorial regions for rapidly rotat- 
ing B stars. Indeed, the rapid rotation leads to polar brighten- 
ing that increases the poles temperature to the hot side of the 
jump. At the same time, the rotation leading to gravity dark- 
ening, the equatorial region may be on the cool side of the 
jump. Consequently, the wind in the equatorial region is denser 
than the wind in the polar region. Nevertheless. lPelupessv et al 



(2000) show that the rotational velocity of the star should be 
very close to its critical value to allow the equatorial wind to 
reach the density needed to create the disk. However, super- 
giant stars cannot be close to critical rotational velocity because 
of probable disruption. Additional mechanisms are therefore 
needed t o allow the sup ergiant star to reach its critical velocity 
(see e.g. lOwockfeOOfjl) . In the particular case of an sgB[e] star 
in an X-ray binary system, the spin-up should occur during the 
supergiant phase of the companion, which indicates a different 
evolutionary stage from other HMXBs. 

This disk itself cannot explain the strong MIR excess the 
sgBTel stars exhibit. Nevertheless, iBjorkman & Cassinelli 
(119931) have shown the existence of a zone in the disk (about 
50-60 stellar radii from the star) in which the temperature 
is below the temperature of sublimation of the dust (about 
1500 K) and the density high enough to allow for its creation. 

IGR J16318-4848 is the source in our sample that exhibits 
the strongest MIR excess, and we believe it is due to the sgB[e] 
nature of its companion star. Indeed, many other strongly ab- 
sorbed sources in our sample do not present any MIR excess. 

Moreover, it is suggested that IGR J16358-4726, the second 
source in our sample that exhibits a MIR excess and whose 
SED needs an additional component to be properly fitted is an 
sgB[e], because its spectrum has all the characteristic features 
of supergiant stars plus the [Fen] feature (CHA08). Our fit is 
therefore in good agreement with their result, and the only other 
source of our sample that definitely exhibits a MIR excess is 
indeed an sgB[e] star. 

In the case of IGR J16195-4945, we are more cautious 
concerning the presence of warm dust that could be responsible 
for a MIR excess, as we lack data above 8 fim. Indeed, Fig 4. 
shows that these source could exhibit a MIR excess, but one 
much lower than the other two. Nevertheless, if this excess 
were to be confirmed, we believe it would be also due to the 
sgB[e] nature of the companion. 

We would like to point out that, because the dust is most lo- 
cated in an equatorial disk in an sgB[e] star, the simple model 
we used to fit the SEDs cannot reproduce the complex distribu- 
tion of the dust around these stars. Nevertheless, it allows the 
detection of a warm MIR excess because of the presence of dust 
in the stellar winds. Finally, for all the stars in our sample, we 
cannot exclude the presence of a cold component - responsible 
for their intrinsic absorption - which we cannot detect because 
of the lack of data above 20 fim. 

4.2. Spectral type and distance 

In our sample, six sources are supergiant stars with a known 
spectral type - IGR J16318-4848 and IGR J16358-4726 are 
sgB[e], IGR J16465-4507 is a B0.5I, IGR J17391-3021 is an 
08Iab(f), IGR J17544-2619 is an 091b, and IGR J19140+0951 
is a B II - and three are found to be O/B supergiants whose 
temperatures derived from our fits allow an asse ssment of the 



spectr al ty pes using the c l assific ation given in iMartins et al 



d2005h and lCrowtheret all J2006) for O and B galactic super- 
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Table 7. Summary of spectral types (SpT) and distances (D») 
derived from our fits for confirmed supergiant stars in our sam- 
ple. (*) sources with an accurate spectral type found in the liter- 
ature (') confirmed supergiant stars whose temperature derived 
from our fits was used to assess their accurate spectral type. 
References to the determination of the spectral type and/or 
spectral class of these sources are found in Table 1 . 



possible to assess the distance of O stars in our sample using 
the standard relation: 



Sources 


SpT 


D,(kpc) 


IGRJ163 18-4848* 


sgB[e] 


~ 1.6 


IGR J 16320-475 1 T 


081 


~ 3.5 


IGR J16358-4726* 


sgB[e] 


~ 18.5 


IGRJ16465-4507* 


B0.5I 


~ 9.4 


IGR J16479-4514 1 " 


08.51 


-4.9 


IGR J17252-3616 r 


08.51 


~ 6.1 


IGRJ17391-3021* 


081 


~ 2.7 


IGR J17544- 2619* 


091 


~ 3.6 


IGR J 19 140+0951* 


B1I 


~ 3.1 



Table 8. Summary of the distances (D„) derived from our fits 
for sources with unconfirmed spectral classes. 



Sources 


SpT 
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4.1 
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giants, respectively, given the uncertainties of observational 
results (~ 2000 K) and uncertainties on the fits temperatures 
as given in Table 6. We therefore found that IGR J16320-4751 
could be an 081 and IGR J16479-4514 and IGR J17252-3616 
could be 08.51 stars. 

Concerning the last three sources whose spectral class is 
unknown, results of the fits listed in Tables 5 and 6 show that 
they are probably all O/B massive stars, and we also used their 
derived temperatures to assess their spectral type using the 
classification given in both papers quoted above. IGR J16195- 
4945 could be a B 1 star, and as already stressed above, it could 
be an sgB[e] due to its MIR excess, IGR J 16207-5 129 and 
IGR J16418-4532 could be 07.5-08.5 stars. Nevertheless, 
even if the high intrinsic X-ray absorption of their associated 
compact objects points towards a supergiant nature since the 
accretion is likely to be wind-fed, the fits themselves do not 
allow an assessment of their spectral classes. We then consider 
they could be either main sequence, giant, or supergiant stars. 



Martins & Pled (120061) give a UB VJHK synthetic photom- 
etry of galactic OI, OIII, and OV stars, with which one can get 
the expected unabsorbed absolute magnitude in J band Mj for 
stars having a given spectral classification. Using the absorbed 
apparent magnitudes mj of our sources and the J band absorp- 
tion we derived from our fits, Aj = 0.289 x Ay, it is then 



D _ 10 0.2(otj-Aj-Mj + 5) 



in pc 



We did not find any synthetic photometry for B supergiants. 
Neverthele ss, expected r adii o f galactic BI, Bill, and BV stars 



are given in lVacca et al.l dl996l) . and we divided these values by 
the jj- ratio derived from our fits to get the star distance. The 
derived distances for the sources whose spectral class is known 
are listed in Table 7, in Table 8 for the others. 



4.3. X-ray properties 

Except in the case of IGR J17544-2619, X-ray absorptions are 
systematically significantly larger than the visible absorptions. 
This indicates the presence in the system of a two-component 
absorbing material: one around the companion star, responsi- 
ble for the visible absorption, and a very dense one around the 
compact object coming from the stellar winds that accrete onto 
the compact object and are responsible for the huge X-ray ab- 
sorption those sources exhibit. 

The obscuration of the compact object by the stellar 
wind is caused by the photoelectric absorption of the X- 
ray emission by the wind, and this absorption varies along 
the orbit of the compact object. This orbital dependence has 
for instance been observed and modelled on 4U 1700-37 
(lHaberl. White. & Kallmanlll989h . 

Moreover, an effect on the X-ray absorption by the 
photoionisation of the stellar wind in the vicinity of the 
compact object by its X -ray emission was predicted by 
Hatchett & McCrayl dl977h . Indeed, in SGXBs, the compact 
object moves through the stellar wind of the companion star, 
and the X-rays are responsible for the enhancement or the 
depletion of the ionised atoms responsible for the acceleration 
of the wind (e.g. Civ and Nv). This has a direct consequence 
on the velocity profile of the wind; when the wind enters into 
an ionised zone, it follows a standard CAK law until it reaches 
a location in which it is enough ionised for no further radiative 
driving to take place, and the wind velocity is "frozen" to a 
constant value from this point. This results in a lower wind 
velocity close to the compact object and consequently a higher 
wind density that leads to a higher obscuration of the compact 
object. 

Most of the sources studied in this work are very absorbed 
in the high-energy domain. Nevertheless, this absorption may 
not be always that high. In the case of very wide eccentric or- 
bits, the column density of the sources could normally vary 
along their orbit and suddenly increase when very close to the 
companion star because of the wind ionisation. In contrast, if 
these objects were to be always very absorbed, it could mean 
that their orbit is very close to the companion star and weakly 
eccentric. If this effect were to be observed, we think it could 
explain the difference in behaviour between obscured SGXBs 
(close quasi-circular orbits) and SFXTs (wide eccentric orbits). 
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Table 9. Sample of parameters we used to fit the SEDs of the 
isolated supergiants. We give their galactic coordinates, their 
spectral types, the interstellar extinction in magnitudes Aj and 
then the parameters themselves: the extinction in the optical 

n 

Ay, the temperature T, and the jy- ratio of the star. 



Sources 


/ 


b 


SpT 


A i 


A v 


T,(K) 


R. 
D. 




HD 144969 
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26000 
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10 
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B0.5Ib 


0.75 
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24700 


8.76 x 10~ 


ii 
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339.38 
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0.81 


1 


24000 


5.30 x 10- 


10 


HD 151804 


343.62 


1.94 


08Iaf 


0.83 


1.3 


32000 


4.20 x 10- 


10 


HD 152234 


343.46 


1.22 


B0.5Ia 


1.17 


1.5 


25100 


4.94 x 10- 


10 


HD 152235 


343.31 


1.1 


Blla 


3.37 


3.9 


24500 


1.13 x 10 


-9 


HD 152249 


343.35 


1.16 


091b 


1.34 


1.7 


30100 


2.89 x 10- 


10 


HD 156201 


351.51 


1.49 


B0.5Ia 


2.68 


2.9 


26500 


2.90 x 10- 
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4.4. Optical properties 

We were able to fit all but three sources with a simple stellar 
black body model. For these three sources, we explained that 
the MIR excess was probably caused by the warm dust created 
within the stellar wind due to the sgB [e] nature of the compan- 
ions. Therefore, it seems that the optical-to-MIR wavelength 
emission of these SGXBs corresponds to the emission of ab- 
sorbed blue supergiants or sgB [e] . 

Moreover, the results of the fits listed in the Table 5 show 
that it is a priori impossible to differentiate an obscured SGXB 
and an SFXT from their optical-to-MIR wavelength SEDs, and 
it then seems that the difference in behaviour between both 
kinds of SGXBs only depends on the geomet ry of the system, 
i.e. it s orbital distance or its orbit eccentricity (IChatv & Rahoui 
2007). 



Nevertheless, to assess a possible effect of the compact ob- 
ject on the companion star, we took a sample of eight isolated 
O/B supergiants in the direction of the Galactic centre and fitted 
their optical-to-NIR wavelength SEDs with an absorbed stel- 
lar black body. The best-fitting parameters are listed in Table 
9 along to their galactic coordinates, their spectral types and 
the interstellar Hi absorption (Aj). The distances of these su- 
pergiants are known, which allowed us to calculate Aj out to 
their position using th e tool available on the MAST website 
dFruscione etal.lll994 . 

We see that their visible absorption is of the same order of 
magnitude as the interstellar H i absorption and well below the 
level of absorption of our sources. This could mean that some 
supergiant stars in SGXBs exhibit an excess of absorption due 
to a local absorbing component. Unfortunately, the total inter- 
stellar absorption out to the distance of our sources is unknown, 
and we cannot compare their visible absorptions derived from 
our fits to the total interstellar absorption out to their position. 

Nevertheless, if this was the case, we think that this ex- 
cess of absorption could also be caused partly by the photoion- 
isation of the wind in the vicinity of the companion star by 
the high-energy emission of the compact object, as this would 



make their winds slower than in isolated supergiant stars. Since 
the wind velocity is lower, the medium is denser and suitable 
for creating a more absorbant material. 

Indeed, in the case of persistent sources with very close and 
quasi-circular orbits, we think that this possible effect could 
be particularly strong, since the wind around the companion 
star would be permanently photoionised and would have lower 
velocities than in isolated supergiants. This could be the general 
scheme of obscured SGXBs. 

On the other hand, in the case of very wide and eccentric 
orbits, the compact object would be most of the time far from 
the secondary and its X-ray emission would not photoionise the 
wind close to the companion star, which would not exhibit any 
visible absorption excess until the compact object got closer. 
This could be the general scheme of SFXTs. 

At last, in both cases, it would be possible to observe a 
variation in the P-Cygni profiles of the companion star (i.e. 
a variation in the wind velocity) with the phase angle of the 
compact object along its orbit. 

As a possible confirmation of this general behaviour, we 
point out that the visible absorptions derived from our fits 
for the companion stars of the only sources in our sample 
that surely exhibit the SFXT behaviour (IGR J16465-4945, 
IGR J17391-3021, and IGR J17544-2619) are far smaller than 
the visible absorptions of the others. Moreover, concerning ob- 
scured SGXBs, the wind velocity of IGR J16318 -4848 was 
found to be ~ 410 km s" 1 dFilliatre & Chatyll2004h . far lower 
than the expected wind velocity for O/B supergiants (~ 1000- 
2000 km s- 1 ). 



5. Conclusions 

In this paper, we presented results of observations performed at 
ESO/VLT with VISIR, which aimed at studying the MIR emis- 
sion of twelve INTEGRAL obscured HMXBs, whose compan- 
ions are confirmed or candidate supergiants. Moreover, using 
the observations performed at ESO/NTT and reported in the 
companion paper (CHA08), previous optical/NIR observations 
found in the literature and archival data from USNO, 2MASS, 
and GLIMPSE, we fitted the broad-band SEDs of these sources 
using a simple two-component black body model in order to 
obtain their visible absorptions and temperatures, and to assess 
the contribution of their enshrouding material in their emission. 

We confirmed that all these sources were likely O/B su- 
pergiant stars and that, for most of them, the enshrouding 
material marginally contributed to the emission. Moreover, in 
the case of IGR J16318-4848, IGR J16358-4726, and perhaps 
IGR J16195-4945, the MIR excess could be explained by the 
sgB[e] nature of the companion stars. 

By comparing the optical and high-energy character- 
istics of these sources, we showed that the distinction 
SFXTs/obscured SGXBs does not seem to exist from optical- 
to-MIR wavelength. Nevertheless, most of the sources in 
our sample are significantly absorbed in the optical, and we 
think that the wind can be denser around some supergiants 
in SGXBs, which could be due to the photoionisation by the 
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high-energy emission of the compact object. 

Several improvements in our study are needed to allow 
definitive conclusions. Indeed, the data used to perform the 
SEDs were not taken simultaneously, which can for instance 
lead to an incorrect assessment of the MIR excess in the 
emission. Moreover, the lack of optical magnitudes for several 
sources could have led to an incorrect fitting of their intrinsic 
visible absorption A v . Finally, the absence of an accurate 
measurement of the total interstellar absorption out to the 
distance of these sources does not allow us to say whether 
the presence of the compact object can lead to a stellar wind 
denser in some supergiants belonging to SGXBs than in 
isolated supergiants. 

We then recommend further optical investigations of these 
sources to study any possible variation in their P-cygni profile 
with the phase of the compact object. We also think that the 
measurement of the distance of these sources is crucial to al- 
low a good assessment of the real interstellar absorption up to 
their distance, in order to detect any local absorbing component 
around companion stars. Finally, we recommend X-ray moni- 
toring so as to study the dependence of their column density 
on orbital phase angle, which could help for understanding the 
difference between obscured SGXBs and SFXTs. 
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Note added i n pro of: After the submission of this pa- 
per, iMoon et aL (2007) reported MIR spectroscopy of IGR 
J16318-4848, confirming the presence of two components of 
warm and cold dust around the sup ergiant B [el companion star. 



Moreover, lHannikainen et alj (|2007) reported a more accurate 



B0.5Ia spectral type of IGR J19140+0951. 



